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Abstract 



Q ■ A class of recently introduced multi-states XX models is generalized to include a defor- 

1^ I mation parameter. This corresponds to an additional nearest-neighbor CC interaction in 

^^ ■ the defining quadratic hamiltonian. Complete integrability of the one-dimensional models 

is shown in the context of the quantum inverse scattering method. The new i?-matrix is 

derived. The diagonalization of the XXC models is carried out using the algebraic Bethe 

Ansatz. 
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1 Introduction 

Multi-states bosonic Hubbard models which generalize the usual su{2) model were recently 
introduced. These models were shown to be integrable and to have an extended symmetry 
H, §, ^ . They are built by coupling two copies of XX 'free-fermions' models ||^, |5|, |3| . These higher 
dimensional 'fermion' models were not obtained as part of a quantum group structure. The 
quantum group approach to the construction of solutions to the Yang-Baxter equation results 
in higher dimensional models with symmetries reflecting the root structure of the underlying 
algebra Q. The combination lowering-raising operators in the i?-matrix generically involves 
all the roots of the algebra. The 'free-fermions' models are also higher dimensional and they 
have simple expressions in terms of the generators of the fundamental representations of su{n). 
However they involve only a subset of the set of roots of su{n). The diagonalization of the 
models reveals an su(2) structure which does not correspond to the known integrable higher- 
spin representations. 

An integrable deformation of the XX models is found and studied in this paper. For the 
lowest dimensional model the XXC model is nothing but the usual su{2) XXZ model. The 
i?-matrix for the higher dimensional models retains the form of the XXZ model and can be seen 
as a multistate generalization of the latter model. The matrix is given and shown to satisfy 
the Yang-Baxter equation. Specific representations are given. Integrability of the models is a 
simple consequence of the quantum inverse scattering framework. The symmetries of the new 
models are obtained. Diagonalization of the XXC models using the algebraic Bethe Ansatz 
method shows that their spectrum is highly degenerate, and is a generalization of the su{2) 
XXZ spectrum. I conclude with some remarks. 

2 The XXC models 

Consider a matrix R of the form 

R[X; j) = p(i) sin 7 + P(2) sin(^ - A) + P(^) sin A (1) 

where 7 is an arbitrary complex parameter and A is the spectral parameter. The operators 
P^^' act on the tensor product of two copies of the same space C". They are taken to satisfy 
the algebra derived in 0]. Let Ei = P^i^i, that is, Ei acts non-trivially on the adjacent spaces 
i, i + 1. The defining relations of the 'free-fermions' algebra A are: 

{Ef,Ei±i} = Ei±i , Ef = E^ (2) 

E,E,±iEi = , EiEj = EjEi for |i - j| > 2 (3) 

where {A, B} = AB + BA. The fourth equation just expresses the fact that Ei and Ej commute 
when they act non-trivially on disjoint spaces. The operators P^^' and P^"^' form a complete set 
of projectors on the tensor product space C" (8) C": 

pW+pC^)=l^ (p{l))2 ^p{l) ^ ^p(2)^)2 ^ p{2) ^ p{l)p(2) ^ p{2)p(l) ^ Q ^4^) 

where I is the identity operator. The operator P^^' is equal to the square of the operator P^^' , 
and p(2)p(3) = p(3)p(2) = Q^ At 7 = 7r/2 the matrix (j^ is equal to the matrix found in g. 
One then checks that R satisfies the Yang-Baxter equation 

RuiX; 7)i?23(A + fi; 7)^i2(^; 7) = ^23(m; 7)-Ri2(A + ^; 7)-R23(A; 7) (5) 

where the parameter 7 is the same for all six factors. The argument 7 of i? will not be written 
in what follows. The calculation is straightforward if tedious. One can expand both sides of the 



YBE, use elementary trigonometric relations and the relations satisfied by the P'-^^'s to conclude. 
Alternatively, setting q = e^^ and y = e and multiplying the YBE by 2iqy, one can also verify 
that the YBE is satisfied for the six values q = 0, ±1, ±i, cx). It is then enough to check that the 
coefficient of q^ vanishes identically to conclude that the YBE is satisfied for all values of A, /i 
and 7. 

It is easy to check that the regularity and unitarity properties hold: 

R{0) = I sin(7) , R{X)R{-X) = I sin(7 + A) sin(7 - A) (6) 

The rational limit of R is obtained by letting A — > 7A, dividing by sin 7 and taking the limit 
7 — > 0. These manipulations conserve all the properties of the ^-matrix. In particular, one 
obtains 

i?(A) = PW + (1 _ A)p(2) + xp(3) (7) 

i?(0)=I, P(A)P(-A) =I(1-A2) (8) 

Explicit representations of the foregoing algebra were given in [||]. Let E""^ be the n x n 
matrix with a one at row a and column /3 and zeros otherwise. Let n, rii and 71-2 be three 
integers such that 

2 < n , 1 < ni < 712 < n — 1 , rii + 112 = n (9) 

and A, B be two disjoint sets whose union is the set of basis states of C", with card(A) = ni 
and card(i?) = 722. Let also 
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The 71-1. ri2 parameters Xa(3 are arbitrary complex numbers. Latin indices belong to A while greek 
indices belong to B. These operators satisfy all the defining relations of the algebra A. The 
restriction ni < ri2 just avoids a double counting of distinct models since one has the obvious 
symmetry ^ <-> P. For rii = 1, ri2 = 72 — 1, and all the Xa/3 equal to each other, one recovers 
the su{n) XX models found in [^. Allowing the twist parameters x to be unequal amounts to 
a multiple deformation of these models. 

It is possible to perform a 'gauge' transformation which puts these models in a Temperley- 
Lieb form. One obtains 

pTL(^) ^ J gij^(^ _ ^) ^ (p(3) ^ pTL) gij^(^) (;L3) 

Py^ = EaeA E/3GB {e'^' E^^ «> E"" + e-'^E"" S^^) (14) 

R^^{X) satisfies the Yang-Baxter equation. 

The transfer matrix is the generating functional of the infinite set of conserved quantities. 
Its construction in the framework of the Quantum Inverse Scattering Method is by now well 
known [0, g, ^. The Lax operator on a chain at site i with inhomogeneity /Xj is given by: 

Loi(A) = Roi{X - fii) = VoiRoiiX - Hi) (15) 



where V is the permutation operator on C" (8) C". The monodromy matrix is a product of Lax 
operators, r(A) = Loa/(A)...Loi('^)) where M is the number of sites on the chain and is the 
auxihary space. The transfer matrix is the trace of the monodromy matrix over the auxihary 
space 0: r(A) = Tro [(Lqa/ ...Lqi) (•^)]- A set of conserved quantities is then given by 

The YBE imphes the following intertwining relations for the elements of the monodromy matrix: 

R{Xi - A2) r(Ai) r(A2) = r(A2) «> r(Ai) r{Xi - A2) (i7) 

Taking the trace over the auxiliary spaces, and using the cyclicity property of the trace, one 
obtains [T(Ai),r(A2)] = 0. The hamiltonians Hp therefore mutually commute. 
The quadratic hamiltonian calculated from ([I^) is equal to 



^^=-^s:(4?.-p.g.c„s,) (18) 

Periodic boundary conditions are implied. These hamiltonians can be written simply in terms 
of su{n) hermitian traceless matrices. For \xa/3\ = 1 and 7 real the hamiltonians are hermitian. 
The rational limit yields: H2 = J^iiPii+i ~ ^u+i)- ^^ ^^ possible to get a positive relative 
sign between P^"^' and P^'^'; one lets 7 — > 7 + vr, and takes the rational limit for the rational 
hamiltonian. 

For all the representations of the operators P*-*^ given above one can introduce the 'conjuga- 
tion' matrix [^] 

C = ^ S^^ - ^ E'"' (19) 

/3g_B aeA 

and verify that p(^) = ^ (I + C (g) C) . Dropping the part proportional to the identity in ( [T^ ) 
gives the quadratic XXC hamiltonian. For n = 2 this is nothing but the XXZ model for su{2). 
For all the parameters equal to one parameter, Xap = x , Va G j4 , V/? G B, the linear 
magnetic-field operators 

Hf = J2Ef a,beA , Hf = Y,Ef a,PeB (20) 

i i 

commute with the transfer matrix. The proof of still holds. The full local symmetry is 
su{ni) © su{n2) 'w(l)- One consequence is that one can add magnetic-field terms for each 
symmetry generators, without spoiling the integrability of the models. 

3 Algebraic Bethe Ansatz 

The diagonalization by algebraic Bethe Ansatz of the XXC models is similar to the one for the 
XX models. We refer the reader to [^, ^ for details and give here the results. I am considering 
the case Xap = x to avoid unessential complications which are peculiar to the higher dimensional 
models (n > 3). 

It is easy to see that the vector ||1) = |1) (8) ... (8) |1) is an eigenvector of the transfer matrix. 
A set of Bethe Ansatz eigenvectors is then obtained from the action of certain elements of the 
monodromy matrix, the C^ = Ti„, and is given by: 

|Ai,...,Ap)= Y. i^"--'"^C„,(Ai)...ajAp) 111) (21) 



where the parameters Aj and the coefficients F are to be determined. This Ansatz also vanishes 
identically iip > M; the proof of this fact is similar to that of [^J. 

One then applies the transfer matrix on the state |Ai, ..., \p) and with the help of relations 
(|l7| ) commutes it through the C^s. The corresponding eigenvalues of t(A) are 

/ M \ p M 

Aiin^,n,),M)^X■,{^,,}) = ijl sm{j - (X - ^,))] ^ f(X^ - X) + l[x-' sm{X - f^,) 

\i=l ) j=l i=l 

P I M \ 

X a("2,p) -Q ;(^ _ ^ .) + (^^ _ 1) Jl ^ sin(A - ^,) bpM (22) 



a 



where /(A) = —-^\-' Here A^"'^:?) jg ^n eigenvalue of the unit-shift operator r'^^'^-*, for 
chain of p sites and 722 possible states at each site. The coefficients i^"p'---'"i are such F is an 
eigenvector of t^^^'^' for the eigenvalue A^"^'^^ Finally the Bethe Ansatz equations are just 

(-ir^ fj -Ml- (A. -/..)) , ^(..„) g sm(7 - (A, - A.)) ^^^^^ ___^^ ^^3^ 

r=i sm(Aj - ^i) ^^. sm(7 - (Afc - Xj)) 

Note that these equations imply the vanishing of the residues of A(("'1'"2)'^^)(A; {/Uj}) at the Aj's. 

The operator r^"^.?) can be written as a product of disjoint permutation cycles. One also 
^(n2,p)\ _ J ^ "Yhe eigenvalues A'^^^'^^ are then q roots of unity, where g is a divisor of p, 
and are highly degenerate. The degeneracy of an eigenvalue depends on p and n2- 

One can perform the above diagonalization procedure over the pseudo-vacuum \\a) (a 7^ 1). 
The set of eigenvalues is exactly the same as the one found above and the eigenvectors have the 
same structure but form a completely distinct set, at least for < p < M. One can also start 
with any of the vectors ||/3), /? € B and obtain yet other sets of eigenvectors. The superscript 
n2 is replaced by ni in (|2^), A'*^^'^^ and r^"-^'^'. These features reflect a large degeneracy of the 
spectrum. 



has (7 



The rational limit is obtained by letting A* -^ jK, l^i -^ Jl^i, dividing the eigenvalue (|2^ ) 
by (sin 7) and taking the limit 7^0. For 7 = 7r/2 the above results reduce to those of |^]. 
For arbitrary 7 the eigenvalues and the Bethe Ansatz equations have the form of those of the 
■sm(2) XXZ model, if one momentarily 'forgets' the A'"'^.?) correction. One can also diagonalize 
the unit-shift operator by algebraic Bethe Ansatz, since this operator appears systematically as 
one of the conserved quantities for all the su{n) XXZ spin-chains. This reflects the underlying 
symmetry seen above. 

A detailed study of the spectrum is desirable. In particular, the thermodynamic limit of the 
spectrum would shed some light on the differences between the XXZ spectrum and the XXC 
one. Taking the logarithm of the Bethe Ansatz equations shows that the distribution of integers 
characterizing the solutions will get a contribution from the A'"'^'^) correction. This may result 
in a denser spectrum. Whether this influences the central charge and the conformal weights 
remains to be clarified. 

A feature of the models with n > 3, is the term 5pM in the expression of the eigenvalues. This 
seems to imply a jump in the spectrum of r(A). However note that, because of the definition 
of the conserved quantities, the hamiltonians Hi, ...,Hm are not affected by this 5-term. For 
vanishing inhomogeneities, the M power of the eigenvalues of r(0)/(sin7) is equal to one 
because one has the unit-shift operator for M sites and n states; one can verify that the M^^ 
power of the eigenvalues (^) is equal to one, after using equations (^) . The eigenvalues of H2 



(at fii = 0) for the trigonometric and rational cases are equal to 

-7^A=o = -M— ^+sm7^-— — 24 

dA sin7 '^ sin(Ajj sin(7 — Ajj 



rflnA,, ,, ^ 1 



|A=0 



"^Iw 



dX '^-^ ^A,(l-A,) 



(25) 



Finally, for 7 7^ tt/2 + ikn, k (z Z, the decorated YBE is not satisfied and it is not possible 
to couple, a la Shastry Q, ^, two XXC copies for generic values of 7. 
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